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In general the experiments on the linear optical properties of a single-layer two-dimensional atomic 
crystal are interpreted by modeling it as a homogeneous slab with an effective thickness. Here I 
fit the most remarkable experiments in graphene optics by using the Fresnel coefficients, fixing 
both the surface susceptibility and the surface conductivity of graphene. It is shown that the 
Fresnel coefficients and the slab model are not equivalent. Experiments indicate that the Fresnel 
coefficients are able to simulate the overall experiments here analyzed, while the slab model fails to 
predict absorption and the phase of the reflected light. 


Graphene, a two-dimensional (2D) carbon-based crys¬ 
tal, was the first atomically thin layer of atoms that was 
produced [i|. The family of 2D crystals however is not 
limited to graphene. Boron nitride or the transition- 
metal dichalcogenides have been grown as atomic mono- 
layers 0. The characteristics of the 2D materials might 
be very different from those of their 3D precursors. This 
is true also for their optical properties. For example, it 
was argued that the high-frequency conductivity (cr) in 
graphene should be a universal constant equal to jAh 
(3|. This is due to the exactly zero overlapping between 
the valence and the conduction bands in graphene (a 
zero-gap semiconductor), while it is finite in graphite. 
Atomically thin transition-metal dichalcogenides are di¬ 
rect band semiconductors while the bulk materials 
have got an indirect band gap. 

Also the optical reflectivity, transmissivity, and ab¬ 
sorption of a 2D crystal are remarkable. The absorption 
of graphene is determined by the fine-structure constant 
0 . Microfabrication of devices used in many experimen¬ 
tal studies currently relies on the fact that 2D crystals 
can be visualized using optical microscopy if prepared 
on top of suitable substrates [Mi- This is due to a 
multiple-reflection effect that enhances the visibility of 
the atomically thin layer, where the optical properties of 
the 2D crystal play a fundamental role. Measurements of 
the optical dielectric functions of 2D crystals by spectro- 
scopic ellipsometry have already been reported m-Ei. 
The interpretation of these results relies on a model of 
the monolayer treated as a homogeneous medium with 
an effective thickness [B-Ej. 

An alternative to this model is to treat the 2D crystal 
as a boundary and to provide the right boundary condi¬ 
tions for a single atomic layer. This approach has been 
used in different papers to derive the reflection and the 
transmission coefficients o f lig ht between two media sepa¬ 
rated by a graphene flake 2ll-l25l|. No reflection or trans¬ 
mission coefficients for other 2D materials different from 
graphene have been considered until now. Although the 
approach used in 21-2^ is substantially correct, these 


papers are primarily interested in the microscopic theo¬ 
retical description of graphene. In general a comparison 
with experimental results of the reflection and transmis¬ 
sion coefficients so obtained was limited to absorption 
and it was not extended to other experiments such as 
optical contrast or ellipsometry. 

Here I deduce an expression for the Fresnel coefficients 
valid for any single-layer 2D atomic crystal. I adopt a 
classical macroscopic approach. I consider first the case 
of a freestanding non-conducting material (such as boron 
nitride 26| ) and then I will turn my attention to optics of 
conducting media (such as graphene). Then will be the 
turn of 2D crystals on substrates. The formulas obtained 
are functional to a comparison with published experi¬ 
mental results on graphene, by far the most studied 2D 
crystal, and they will enlighten the power of Fresnel co¬ 
efficients for simulating the overall phenomena described 
above. 

Consider a flat 2D crystal, composed of N atoms per 
square cm with an atomic polarizability a 27| . If we ap¬ 
ply an electric field in the plane of the crystal a macro¬ 
scopic dipole moment arises and it is possible to define a 
density of polarization P. If the electric field is applied 
orthogonally to the 2D crystal, no macroscopic polariza¬ 
tion is created (or in any case it is much smaller and I do 
not consider it). Indeed to have a macroscopic polariza¬ 
tion the microscopic dipoles need to be aligned, to gen¬ 
erate a macroscopic separation of charges. As a further 
simplification I suppose that the 2D crystal is isotropic 
in its own plane. This seems realistic because graphite, 
for instance, is a uniaxial crystal with the optical axis 
along the graphene’s exfoliation axis. I assume also that 
P = eoxE where eo is the vacuum permittivity and x 
is the electric susceptibility. Wherever the polarization 
in matter changes with time there is an electric current 
Jp, a genuine motion of charges. The connection be¬ 
tween rate of change of polarization and current density 
is Jp = dP/dt. It is important to note that in passing 
from a bulk to a 2D crystal the dimensions of P pass 
from dipole moment/volume to dipole moment/area, the 
dimensions of x pass from a pure number to a length 
and Jp is a current per unit length. It is not the aim of 
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FIG. 1. Oblique incidence of a plane wave on graphene. The 
electric (magnetic) field for s (p) polarization is shown. Inset: 
Three layer substrate; graphene is deposited at the interface 
of media 1 and 2. 


this paper to give a microscopic theory that furnishes y. 
Experiments will fix it. 

Suppose that a 2D crystal is suspended in vacuum and 
that a plane wave falls onto it (Fig. 1; m = n 2 = 1); the 
relation in between E and H in the incident, reflected 
and transmitted waves is r]H = s A E where s is the 
unit vector in along the propagation direction and p is 
the impedance of vacuum. The boundary conditions are 
k A (E 2 — El) = 0, K A (H 2 — Hi) = Jp where k is the 
unit vector along the 0 axis. For s (p) polarization then, 


^xi Exr 
^xi Ext 


— Ext t 
= 

eox’ 


Eyi Hyx - Eyt ~\~ iU^Px'l 


{Eyi Eyr — Eyt) (1) 

{Eyi — Eyr = —-^) 

(^Exi Ext — Ext iujPy^ 


where time dependence exp(ia;t) {ui is the angular fre¬ 
quency of the light) has been assumed. 

Defining 0 = Er/Ei, tg = Et/Ei and Vp = Hr/Ei, 

tp = E[t/Ei as the reflection and the transmission coeffi¬ 
cients, we have 

^ -ikx ^ ikx cos Oj 

ikx+ ‘2 cos 0i’ ikxcos9i + 2 

and ts = rg + 1, tp = 1 — Tp {k is the wave vector of 
light in vacuum and 6i is the angle of incidence). From 
energy flux considerations [ 2 ^, in this special case, the 
reflectivity is = ks(p)li the transmissivity is T's(p) = 
|ts(p)l, their sum Rg[p) -I- T)i(p) = 1 shows that there 
is no absorption. 

We turn now our attention to conducting media. The 
Ohm’s law J = uE must be taken into account; again we 
assume that J can not exist (or it is small) in a direction 
orthogonal to the crystal plane, and in-plane isotropy. 
The boundary conditions for H changes into k A (H 2 — 
Hi) = Jp -I- J, and we add to Eqs. (1) the Ohm’s law: 


Ex 


{Eyi Eyr — 


( 3 ) 


for s (p) polarization. We obtain 

ikx Ear] {ikx + uni) cos 9 1 , 

tp — __ * p (^) 

* ikx + err] + 2 cos 9i' ^ {ikx + enj) cos 9i + 2 

and tg = rg+l^tp = l—Vp. As for a bulk material conduc¬ 
tivity is connected with the transformation of part of the 
electromagnetic energy into heat. At normal incidence, 

4ctp 

Tg(p) + Rgip) = 1 - 4 + 4 ^^ + ^ 2^2 + - 

In the case of graphene, at wavelengths A where a = 
/Ah, we retrieve the remarkable result that the fine- 
structure constant defines its optical transparency Q. 

The reflection and the transmission coefficient for the 
case of a 2D crystal at the interface of two different bulk 
media (1 and 2 in Fig. 1) is now easily obtained. Only 
the relation in between E and H in the incident, reflected, 
and transmitted waves changes: 

Ei(r) ^i(r) A Ei(r\^ St A Ef^ (5) 

ni ' ^ ^ ^ ' 712 

where ni, 712 are the refractive indexes. We obtain: 


Til cos 9i — 712 cos 9t — ikx — urj 

Til cos 9i + 712 cos 9t + ikx + urj 

712 cos 9i — Til cos 9t + {ikx + ut]) cos 9i cos 9t 

712 cos 9i + m cos 9t + {ikx + urj) cos 9i cos 9t 


and tg = Cg + 1, tp = (1 — rp)n 2 cos9i/{ni cos9t). 

The reflection coefficient of a 2D crystal on a strati¬ 
fied medium is of primary importance also. I consider 
the case represented in Fig. 1 (inset), corresponding to 
published experimental work that will be considered in 
this article. The boundary conditions for the interface 
in between bulk media 2 and 3 are k A (E 3 — E 2 ) = 0, 
k A (H 3 — H 2 ) = 0. Equation (5) is easily extended to 
medium 3 [ 2 ^. For s polarization. 


^xi “1" ^xr — ^x-\- ^x—1 
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^xi “1“ ^xr — ; Hyj- — ^y-\- ^y- 

Ex+e-^^ + Ex-P^ = Ext; 
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where P = kn2hcos92, 6*2 is the propagation angle in 
medium 2, and h is its thickness [29|. This furnishes 

^ D2 s + r23.e-^^^ - ti2«fg^(l -f r23se-^^^) 

l + ri2«r23.e-2*/3+ti2,f^(l + r23«e-2*/3) 
ri 2 p + r2Spe-^^P + ti 2 p - r2Spe-^^P) 

” 1 + ri2pr23pe-2*/5 + ti^p - r2Spe-^^P) 
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FIG. 2. The first three graphs show the optical contrast of 
graphene on Si02/Si. The experimental data and the simu¬ 
lations based on the slab model are extracted from [^. The 
last graph shows the absorption of a free-standing graphene 
layer. Experimental data are extracted from 


If X Slid a are zero we obtain the usual reflection coeffi¬ 
cients for a three-layered medium [^. The sign differ¬ 
ence with respect to [2^ in the face of 2i/3 comes from 
the different choice of the temporal dependence 

We can now compare some of the most remarkable ex¬ 
periments in graphene optics with the theoretical predic¬ 
tions given by the Fresnel formulas (4), (6), (9) and with 
the predictions furnished by the slab-based model. All 
the experimental data reproduced in this paper have been 
extracted from the original papers via software digitiza¬ 
tion, ensuring a very precise reproducibility. Throughout 
this paper data are represented by dots, slab model pre¬ 
dictions by a dashed line, and Fresnel model predictions 
by a solid line. The error bars reported in the graphs 
have a length of 2 standard deviations s (see Table I) for 
the respective sets of experimental data. 

As explained before the visibility of graphene is en¬ 
hanced if prepared on top of Si02/Si wafers. The con¬ 
trast is defined as the relative intensity of reflected light 
in the presence and absence of graphene Figure 2 
compares the experimental data and the slab model pre- 
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FIG. 3. Spectral dependance of the ellipsometric parameter A 
for graphene on an amorphous quartz substrate. Experimen¬ 
tal data for the four 9i are extracted from El- The Fresnel 
and the slab model simulations for 9i =60° are reported (see 
text). 


dictions published in @ with the Fresnel theory. The 
first three graphs give the optical contrast for single-layer 
graphene on top of Si 02 /Si wafers with three different 
Si02 thicknesses at 9i = 0. For the optical constants 
of Si02 I used the Sellmeier equations, for Si I made 
a fit from data in (^ . Polarizations s and p give the 
same results. The Fresnel fit reported in the first three 
graphs of Fig. 2 is for cr = e^/4?i = 6.08 ■ 10“® Q, and 
X = 5 ■ m. This is compared with the slab fit used 

in 0 . Anyway the contrast measurements are not able to 
discriminate very well the value of x- All the Fresnel fits 
with a = e^/4?i and x ^ 5 • 10“^° m give in practice the 
same result. We will see that ellipsometric measurements 
solve this problem. 


To quantify how much better the Fresnel-based model 
is, compared to the slab-based one, a chi-squared test 
is reported in Table I where Xof(S) reduced chi- 


squared value for the Fresnel (slab) fit and Pf{S) is the 
probability of getting a value as large as Xof(S) assuming 
that the Fresnel (slab)-based model is correct 3l|. For 
the contrast measurements the standard deviation s has 
been chosen as the diameter of the experimental dots 
in Fig. 2 of Ref. From Fig. 2 and Table I the 

Fresnel-based model is clearly superior to the slab-based 
one. It seems anyway that the predictions of the spectral 
positions of the maxima for h = 300 nm and h = 90 nm in 
Fig. 2 are better guessed by the slab-based model. This 
is also the reason why Xqf i® good for h = 300 nm. 
From simulations it turns out that the positions of the 
maxima are more sensitive to the substrate parameters 
than to the graphene layer; this is a possible explanation 
of the discrepancy from the Fresnel-based ht. 
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FIG. 4. Spectral dependance of the ellipsometric parameter A 
for graphene on an a morp hous quartz substrate for each of the 
four 8i considered in [2. Experimental data are represented 
by dots, slab model predictions by a dashed line and Fresnel 
model predictions by a solid line. 



FIG. 5. Spectral dependance of ellipsometric parameter 4' for 
graphene on an amorphous quartz substrate. Experimental 
data for the four di are extracted from 0. The Fresnel 
and the slab model simulation for each 6i and for A > 350 
nm are reported. Fits are obviously relative to their nearest 
experimental data points. 


The last graph of Fig. 2 considers a free standing 
graphene layer and compares Rs{p) +Tf^) for the present 
theory and for the slab model used in [g . The two theo¬ 
retical predictions are very different. For a constant value 
of (T = /Ah the Fresnel theory predicts a constant ab¬ 
sorption as a function of the wavelength whereas the slab 
model predicts a wavelength dependance. The two the¬ 
oretical predictions are compared with the experimental 
data published in Q. In this case it was possible to re¬ 
trieve s (Table I) from the spreading of the experimental 
data. Again the superiority of the Fresnel model is quite 
evident. 

To fix Xi I turned my attention to spectroscopic ellip- 
sometry of graphene flakes located on a flat amorphous 
quartz. Figures 3, 4, and 5 compare the theoretical pre¬ 
dictions for the Fresnel-based model and for the slab- 
based model with the measurements published in [l5| . 
Figures 3 and 4 show the simulated spectral dependence 
of the ellipsometric parameter A [3^ at four angles of 
incidence. Figure 3 shows the ensemble of the experi¬ 
mental data and the fits of the two models for 9i =60°. 
For clarity reasons Fig. 4 focuses on each 9i. From Figs. 
3 and 4, for A > 350 nm, a good Fresnel fit is obtained for 
a = /Ah and y = 1.0 • 10“® m. This is compared with 
the slab-based model used in 0 in the same frequency 
range. For 9i =55° A does not fit very well, maybe be¬ 
cause it is too close to Brewster and cross polarization 
effects [0 may be present (10“^ < Rp < 1.5 • 10“^ at 9i 
=55° in this spectral range). From Fig. 4 and Table 1 
the Fresnel-based model fits much better than the slab- 
based model the experimental data for the ellipsometric 
parameter A, and hence it better predicts the phase of 


the reflected light. As for the measurements of contrast 
the standard deviations s for the ellipsometric measure¬ 
ments have been extracted from the experimental line 
widths in Fig. 5(b) of Ref. [H. 

The ellipsometric parameter A is very sensitive to the 
graphene film because A = 180 ° or 0° for the quartz 
substrate and all of the non trivial phase contribution 
to the reflection coefficients comes from graphene. In 
particular from formulas ( 6 ) the phase of the reflected 
light is different from 0° or 180° only if y 7 ^ 0 m. The 
experimental data for A allow one to fix a value of x 
within ±1 • 10 “^° m, and they are also sensitive to its 
sign. 

Figure 5 shows the ellipsometric parameter Ik ( 3 ^ . 
From Fig. 5 and Table I the Fresnel fit and the slab 
fit for Ik give almost equivalent predictions, with a slight 
superiority of the slab-based one. In this case the role of 
the dielectric substrate is overwhelming (more than 98% 
of the 41 signal is due to it). This is a possible expla¬ 
nation of the small discrepancy of the experimental data 
from the Fresnel-based fit. The fit reported in Fig.5 is 
for the same values used to fit A. A better fit is obtained 
if we let a be 30% or 40% less than this value. This does 
not really affect A but it can be an indication of defects 
in the sample. In fig 5 dots have the linear dimensions of 
one standard deviation s (experimental line width in Fig. 
5(a) of Ref. 0), this is why error bars are not reported. 


For A < 350 nm data indicate that a and x vary sen¬ 
sibly with A. In this spectral range it is possible for each 
wavelength to fix a different value of cr and x- We con¬ 
sider for instance A = 270 nm. At this wavelength by 
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TABLE I. Chi-squared test value for the Fresnel-based model 
fit and the slab-based model fit; n: number of experimental 
data points considered for the fitting procedure. 




n 

s 

XOF 

Pf 

Xqs 

Ps 

Gontrast 300 nm 

12 

0.01 

4.35 

< 0.1% 

4.93 

< 0.1% 

Contrast 200 nm 

12 

0.01 

0.15 

99.8% 

1.68 

9.7% 

Contrast 90 nm 

12 

0.01 

1.18 

30.2% 

6.42 

< 0.1% 

Absorption 

34 0.0015 

1.09 

33.4% 

6.75 

< 0.1% 

A 

60° 

25 

0.8° 

1.60 

5% 

45.74 

< 0.1% 

A 

55° 

20 

0.8° 

59.87 

< 0.1% 333.37 

< 0.1% 

A 

50° 

23 

0.8° 

0.82 

69% 

18.42 

< 0.1% 

A 

45° 

28 

0.8° 

0.55 

94% 

7.74 

< 0.1% 

4- 

60° 

29 

0.4° 

0.98 

49% 

1.92 

0.4% 

T 

55° 

29 

0.4° 

2.21 

< 0.1% 

2.48 

< 0.1% 

4' 

50° 

29 

0.4° 

0.96 

53% 

0.77 

79% 

4' 

45° 

29 

0.4° 

0.78 

77% 

0.47 

98% 


setting a = 3.1 • /Ah and x = —1-2 -10 ® m we ob¬ 
tain for 6i = 60°, 55°, 50°, 45° respectively = 3.2°, 
4.8°, 12.1°, 18.9° and A = 18.2°, 168°, 175°, 177.5°, and 
an absorption of 6.5%, in very nice agreement with data 


sorptioi 

in Ref. [15|. From simulations: |x| is still too small to 
influence absorption and so it is fixed by ct, and x fixes 
the rest; A is sensitive to the sign of X) 'I' to its magni¬ 
tude only. A negative value of x at 270 nm is probably 
related to the fact that the real part of the dielectric 
constant of graphite become negative around 4 eV and 
up to around 7 eV [3^. Fresnel coefficients simulate as 
well spectroscopic ellipsometry of graphene deposited on 
Si02/Si reported in the same paper 1^. I confirm that 
4/ is sensitive to presence of the graphene layer. Unfor¬ 
tunately (even up to |x| = 5 • 10“® m) A is dominated by 
the signal from the substrate. 


In conclusion the comparison with the experimental 
results shows that the Fresnel coefficients are essential to 
interpret the most remarkable experiments in graphene 
optics. Any hypothesis on an effective thickness of a 
single-layer 2D atomic crystal as required by modeling it 
as an homogeneous slab is not necessary. This last model 
is not able to reproduce properly either the absorption 
of graphene, or the phase of its reflection coefficient. On 
this basis any physical parameter, deduced from it, is 
hardly meaningful. As for bulk materials, ellipsometry 
is able to furnish both x and cr, showing that these are 
the physically meaningful quantities experimentally ac¬ 
cessible from the linear optical response of a 2D atomic 
crystal. In particular for graphene, from the ensemble of 
the experimental data considered, in the spectral range 
450 nm < A < 750 nm, cr = 6.08 ■ lO”® ± 2 ■ lO-^D-^ and 
X = 8 • 10“^° ± 3 • 10“^° m. Of course the Fresnel coeffi¬ 
cients are not valid for 2D crystal bilayers or multilayers. 
In these cases light propagates from one layer of atoms 
to the other and a thickness should be considered. This 
shows once more how special single-layer atomically thin 
2D crystals are. 
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